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ABSTRACT 

The hadroproduction experiments HARP, MIPP and NA61 (SHINE) as well 
. as their implications for neutrino physics are discussed. HARP measurements 

^vq ' have already been used for precise predictions of neutrino beams in K2K and 

£H ■ MiniBooNE/SciBooNE experiments and are also being used to improve the at- 

mospheric neutrino flux predictions and to help in the optimization of neutrino 
H" ~i . factory and super-beam designs. The MIPP experiment has nearly finalized 

' measurements of hadron production from the NuMI target used in the MINOS 

(N| . experiment. First measurements released recently by the NA61 (SHINE) ex- 

periment are of significant importance for a precise prediction of the J-PARC 
^ . neutrino beam used for the first stage of the T2K experiment. All three exper- 

QJ ' iments - HARP, MIPP and NA61 - provide also a large amount of input for 

■ validation and tuning of hadron production models in Monte-Carlo generators. 

1. The HARP experiment 

The HARP experiment ID at the CERN PS was designed to make measurements 
of hadron yields from a large range of nuclear targets from hydrogen to lead and 
for incident particle momenta from 1.5 GeV/c to 15 GeV/c. The main motivations 
■^J- ■ are the measurement of pion yields for a quantitative design of the proton driver 

\Q ■ of a future neutrino factory^, a substantial improvement in the calculation of the 

atmospheric neutrino flux ^' and the measurement of particle yi elds as input for the 
flux calculation in accelerator neutrino experiments, such as K2K^^, MiniBooNE^ 1 
^ \ and SciBooNE^ 1 . In addition to these specific aims, the data provided by HARP are 

valuable for validation and tuning of hadron production models used in simulation 
programs. The HARP experiment is described in more details below in order to 
illustrate some common features of modern hadron production experiments. 

To provide a large angular and momentum coverage of the produced charged par- 
ticles the HARP experiment makes use of a large-acceptance spectrometer consisting 
of a forward and large-angle detection system. A detailed description of the exper- 
imental apparatus can be found in Ref. The forward spectrometer — based on 
large area drift chambers ® and a dipole magnet complemented by a set of detectors 
for particle identification (PID): a time-of-flight wall^ (TOFW), a large Cherenkov 
detector (CHE) and an electromagnetic calorimeter - - covers polar angles up to 
250 mrad which is well matched to the angular range of interest for the measurement 
of hadron production to calculate the properties of conventional neutrino beams. The 
large-angle spectrometer — based on a Time Projection Chamber (TPC) located in- 
side a solenoidal magnet — has a large acceptance in the momentum and angular 



range for the pions relevant to the production of the muons in a neutrino factory. 
It covers a large majority of the pions accepted in the focusing system of a typical 
design. The neutrino beam of a neutrino factory originates from the decay of muons 
which are in turn the decay products of pions. 

A full set of data collected by the HARP experiment with thin (5% of nu- 
clear interaction length, Ai) and cryog enic targets have been analyzed and pub- 
lished i m | l^U3 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22t Thoge regultg coyer &u thg physicg sub j ects 

discussed above. 



1.1. Results obtained with the HARP forward spectrometer 

A detailed description of established experimental techn iques for th e data analysis 
in the HARP forward spectrometer can be found in Refs. I10 | H | 12 [I3^ Only a brief 
summary is given here. 

The absolutely normalized double-differential cross-section for a process like p + 
Target — > n + +X can be expressed in bins of pion kinematic variables in the laboratory 
frame (momentum, p n , and polar angle, 9 n ), as 

where is the cross-section in cm 2 /(GeV/c)/srad for each {p^^O^) bin covered 

in the analysis; j^- is the reciprocal of the number density of target nuclei; t is the 
thickness of the target along the beam direction; Ap and AQ are the bin sizes in 
momentum and solid angle (Ap = p max - p min ; AQ = 2ir(cos(6 min ) - cos(9 max )))\ 
iVpot is the number of protons on target after event selection cuts; N n (p n , 6^) is the 
yield of positive pions in bins of true momentum and angle in the laboratory frame. 
Eq. ([T]) can be generalized to give the inclusive cross-section for a particle of type a 

where reconstructed quantities are marked with a prime and M~} , n , / is the inverse of 

1 1 pOap a 

a matrix which fully describes the migrations between bins of true and reconstructed 
quantities, namely: lab frame momentum, p, lab frame angle, 9, and particle type, a. 

There is a background associated with beam protons interacting in materials other 
than the nuclear target (parts of the detector, air, etc.). These events are subtracted 
by using data collected without the nuclear target in place after normalization to the 
same number of protons on target. This procedure is referred to as the 'empty target 
subtraction'. 

The event selection is performed in the following way: a good event is required to 
have a single, well reconstructed and identified beam particle impinging on the nu- 
clear target. A downstream trigger in the forward trigger plane (FTP) is also required 



to record the event, necessitating an additional set of unbiased, pre-scaled triggers 
for absolute normalization of the cross-section. These pre-scaled triggers (e.g. 1/64) 
are subject to exactly the same selection criteria for a 'good' beam particle as the 
event triggers allowing the efficiencies of the selection to cancel, thus adding no addi- 
tional systematic uncertainty to the absolute normalization of the result. Secondary 
track selection criteria have been optimized to ensure the quality of the momentum 
reconstruction as well as a clean time-of-flight meas ureme nt while maintaining high 
reconstruction and particle identification efficiencies I H I -^t _ 

The first HARP physics publication reported measurements of the tc + produc- 
tion cross-section from a thin 5% Ai aluminum target at 12.9 GeV/c proton momen- 
tum. This corresponds to the energies of the KEK PS and the target material used 
by the K2K experiment. The results obtained in Ref. ^® were subsequently applied 
to the final neutrino oscillation analysis of K2K — , allowing a significant reduction of 
the dominant systematic error associated with the calculation of the so-called far-to- 
near ratio from 5.1% to 2.9% (see and ^» for a detailed discussion) and thus an 
increased K2K sensitivity to the oscillation signal. 

The next HARP goal was to contribute to the understanding of the MiniBooNE 
and SciBooNE neutrino fluxes. They are both produced by the Booster Neutrino 
Beam at Fermilab which originates from protons accelerated to 8.9 GeV/c by the 
booster before being collided against a beryllium target. A fundamental input for the 
calculation of the resulting flux is the measurement of the tt + cross-sections from 
a thin 5% Ai beryllium target at exactly 8.9 GeV/c proton momentum. 

The double-differential cross-section for the production of positively charged pi- 
ons from collisions of 8.9 GeV/c protons with beryllium have been measured in the 
kinematic range from 0.75 GeV/c < p n < 6.5 GeV/c and 0.030 rad < 9 n < 0.210 rad, 
subdivided into 13 momentum and 6 angular bins^l. A typical total uncertainty 
of 9.8% on the double-differential cross-section values and a 4.9% uncertainty on the 
total integrated cross-section are obtained. These HARP results have been used for 
neutrino flux predictions in the MiniBooNE and SciBooNE experiments^ 1 . 

Sanford and Wang have developed an empirical parametrization for describing 
the production cross-sections of mesons in proton-nucleus interactions (e.g. p+A — > 
7r + + X). This SW parametrization has the functional form: 

d?a p c ^ p 

:(P, 0) = ex P[ c i - c 3^i c 6 6{p - c 7 p beam cos C8 6)}p C2 {l ) , (3) 



where X denotes any system of other particles in the final state; Pbeam i s the pro- 
ton beam momentum in GeV/c; p and 9 are the ir + momentum and angle in units 
of GeV/c and radians, respectively; d 2 a /(dpdQ) is expressed in units of mb/(GeV/c 
sr); dQ = 2ir d(cos8); and the parameters Ci, . . . , c% are obtained from fits to meson 
production data. HARP data have been fitted using this empirical SW parametriza- 
tion. This is a useful tool to compare and combine different data sets. In a 



harp p-c -> tt^- HARP p-N2 -> tt + - HARP p-02 -> jt + ~ 




Figure 1: Measurements of the double-differential production cross-sections of tt + (open circles) 
and 7T~ (closed circles) from 12 GcV/c protons on carbon (left), N2 (middle) and O2 (right) target 
as a function of pion momentum, p, in bins of pion angle, 9, in the laboratory frame. Different 
panels show different angular bins. The error bars shown include statistical errors and all (diagonal) 
systematic errors. The curves show the SW parametrization of Eq. (JU) with parameters given in 
Ref. 13 - ) (solid line for tt + and dashed line for 



global SW parametrization is provided for forward production of charged pions as 
an approximation of all the studied datasets. It can serve as a tool for quick yields 
estimates. 

The next HARP analysis was devoted to the measurement of the double-differential 
production cross-section of tt^ in the collision of 12 GeV/c protons (see Fig. [1]) and 
pions with a thin 5% Ai carbon target These measurements are important for a 
precise calculation of the atmospheric neutrino flux and for a prediction of the devel- 
opment of extended air showers. Simulations predict that collisions of protons with a 
carbon target are very similar to proton interactions with air. This hypothesis could 
be directly tested with the HARP data. Measurements with cryogenic O2 and N2 
targets^*, also shown in Fig. [Q confirm that p-C data can indeed be used to predict 
pion production in p-0 2 and p-N 2 interactions. 

It is important to emphasize that a systematic campaign of measurements has been 
performed by HARP: all thin target data taken with pion beams are analyzed and 
published now^^. Similar results have also been obtained for incoming protons ^D. 
Moreover, results on forward proton production in p-A and ir^-A interactions have 
been published recently 1 ^. Fig. |2] shows the dependence on the atomic number A of 
the forward proton yields in p-A and ir^-A interactions averaged over two forward 
angular regions (0.05 rad < 9 < 0.15 rad and 0.15 rad < 6 < 0.25 rad) and four 
momentum regions (0.5 GeV/c < p < 1.5 GeV/c, 1.5 GeV/c < p < 2.5 GeV/c, 
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Figure 2: The dependence on the atomic number A of the forward proton production yields in 
p-A and ir^-A (A = Be, C, Al, Cu, Sn, Ta, Pb) interactions averaged over two forward angular 
regions (0.05 rad < 9 < 0.15 rad and 0.15 rad < 8 < 0.25 rad) and four momentum regions 
(0.5 GeV/c < p < 1.5 GeV/c, 1.5 GcV/c < p < 2.5 GeV/c, 2.5 GeV/c < p < 3.5 GcV/c and 
3.5 GeV/c < p < 4.5 GeV/c), for the four different incoming beam momenta (3, 5, 8 and 12 GeV/c). 



2.5 GeV/c < p < 3.5 GeV/c and 3.5 GeV/c < p < 4.5 GeV/c), for the four different 
incoming beam momenta (3, 5, 8 and 12 GeV/c). 

The advantage of all these measurements is that they are performed with the same 
detector, thus related systematic uncertainties are minimized. 

1.2. Results obtained with the HARP large-angle spectrometer 

First results on the measurements of the double-differential cross-section for the 
production of charged pions in p-A collisions emitted at large angles from the incoming 
beam direction have been obtained for the full set of available nuc lear targ ets using an 
initial part of the accelerator spill to avoid distortions in the TPC^ ^^^ . The pions 
were produced by proton beams in a momentum range from 3 GeV/c to 12 GeV/c 
hitting a target with a thickness of 5% Ai. The angular and momentum range covered 
by the experiment (100 MeV/c < p < 800 MeV/c and 0.35 rad < 9 < 2.15 rad) is 
of particular importance for the design of a neutrino factory. Track recognition, 
momentum determination and particle identification were all performed based on the 
measurements made with the TPC (see^^ for more details). 

After the development of a special algorithm to correct for the TPC dynamic 
distortions ^SJ an d validation of the TPC performance with benchmarks based on real 
data ^3, full spill data have been analyzed and published for incoming protons 



HARP 



HARP 



♦ 


pPb 7T + 










1000 





pPb TT — 








-- 


pTa 7T + 










: □ 


pTa 7r~ 










pS n 7T + 












' V 


pSn 7T~ 










pCu 7T + 








O 






pCu 7T~ 
















ft 

▼ 


GeV, 










ft 


ft 


- • 


pC 7T + 




* 




- 


pC TT~ 










pBe 7r + 




▼ 




S-. 






pBe tt~ 


ft 








* 








100 














▼ 






























rO 








/IS 




xx 




▼ 




X 












X 


X 




A. 


X 






































CD 


CD 




• 


* * 


• 


C2-, 








9 


* * 






• 


* 






























fc> 


10 
















O- 

o 


IS p <0. r ?' 

35a e <i 


(GoV/c) 
55 (rad) 


-a 








0. 



IS p <O.T 

35a e <i 


(GcV/c) 
55 (rad) 



Pbe 



5 IO 15 

m [GeV/c] 



5 IO 15 

m [GeV/c] 



Figure 3: The dependence on the beam momentum of the ir + (left) and ir~ (right) production 
yields in p— A (A = Be, C, Al, Cu, Sn, Ta, Pb) interactions averaged over the forward angular region 
(0.350 rad < 8 < 1.550 rad) and momentum region 100 MeV/c < p < 700 MeV/c. The results are 
given in arbitrary units, with a consistent scale between the left and right panel. Data points for 
different target nuclei and equal momenta are slightly shifted horizontally with respect to each other 
to increase the visibility. 



and pions ^* . Fig. [3] illustrates the dependences of the measured charged pion yields 
on the beam momentum in p-A (A = Be, C, Al, Cu, Sn, Ta, Pb) interactions averaged 
over the forward angular region (0.350 rad < 6 < 1.550 rad) and momentum region 
100 MeV/c < p < 700 MeV/c. 

An additional analysis has also been performed on the comparison of the pion 
yields measured with the short (5% Ai) and long (100% Ai) nuclear targets^ 1 . 

All HARP measurements described above have been compared with predictions 
of Monte-Carlo models available within CK.W'T l - !,: and MARS : * (,: . as well as with 
the GiBUU mod el ^D. Some models provide a reasonable description of HARP data 
in some regions I 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 }^ w hji e there is no model which would 
describe all aspects of the data. 

Note that tables with HARP results in text format are available from the DURHAM 
database 123. 



2. The MIPP experiment 

The MIPP experime nt at F ermilab took data in 2005. Some of these data are now 
analysed and published I 33 | 34 | rp^g ana ]y S j s Q f hadron yields in 120 GeV/c proton 
interactions with the NuMI carbon target is nearly complete. These results could be 
used to improve precision of neutrino flux predictions in the MINOS experiment 



There is also a recent proposal for an upgrade of the MIPP experiment 



3. The NA61 (SHINE) experiment 

The physics program of the NA61 or SHINE (where SHINE = SPS Heavy Ion 
and Neutrino Experiment) experiment at the CERN SPS consists of three main 
subjects. In the first stage of data taking measurements of hadron production in 
hadron-nucleus interactions needed for neutrino (T2K ^H) and cosmic-ray (Pierre 
Auger ^3 an d KASCADE^3) experiments are performed. At later stages of the 
NA61 experiment hadron production in proton-proton and proton-nucleus interac- 
tions needed as reference data for a better understanding of nucleus-nucleus reactions 
will be studied, energy dependence of hadron production properties will be measured 
in p-p and p-Pb interactions as well as in nucleus-nucleus collisions, with the aim to 
identify the properties of the onset of deconfinement and find evidence for the critical 
point of strongly interacting matter. 

The NA61 (SHINE) apparatus is a large acceptance spectrometer at the CERN 
SPS for the study of the hadronic final states produced in interactions of various 
beam particles (71*, p, ions) with a variety of fixed targets at the SPS energies. The 
main components of the current detector were constructed and used by the NA49 
experiment . The main tracking devices are four large volume TPCs. Two of them, 
the vertex TPCs (VTPC-1 and VTPC-2), are located in the magnetic field of two 
super-conducting dipole magnets (maximum bending power of 9 Tm) and two others 
(MTPC-L and MTPC-R) are positioned downstream of the magnets symmetrically 
with respect to the beam line. One additional small TPC, the so-called gap TPC 
(GTPC), is installed on the beam axis between the vertex TPCs. The setup is 
supplemented by time of flight detector arrays two of which (ToF-L/R) were inherited 
from NA49 and can provide a time measurement resolution of a to f < 90 ps. A 
new forward time of flight detector (ToF-F) was constructed in order to extend the 
acceptance of the NA61 (SHINE) setup for pion and kaon identification as required 
for the T2K measurements. The particle identification in NA61 is based on the 
differential energy loss dE/dx measured in the TPCs combined with the mass squared 
measurements based on the ToF information, see Fig. HI 

Two carbon (isotropic graphite) targets were used for T2K-related measurements: 
1) a 2 cm-long target (about 4% Ai) with density p = 1.84 g/cm 3 , the so-called thin 
target; 2) a 90 cm long cylinder of 2.6 cm diameter (about 1.9 Ai), the so-called T2K 
replica target with density p = 1.83 g/cm 3 . 

3.1. First results 

The NA61 (SHINE) experiment was approved at CERN in June 2007. The first 
pilot run was performed during October 2007. In total about 670 k events with 
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Figure 4: Examples of PID capabilities of the NA61/SHINE spectrometer for positively charged 
particles. (Left) Specific energy loss in the TPCs as a function of momentum. Curves show pa- 
rameterizations of the mean dE/dx calculated for different particle species. (Middle) Mass squared, 
derived from the ToF-F measurement and the fitted path length and momentum, versus momen- 
tum. The lines show the expected mass squared values for different particles. (Right) Example 
of two-dimensional m 2 -dE / dx plot for the momentum range 2-3 GeV/c. Four clear accumulations 
corresponding to positrons, pions, kaons and protons are observed. 2a contour around fitted pion 
peak is shown. The combination of both m 2 and dE/dx measurements provides close to 100% purity 
in the pion selection over the whole momentum range. 



the thin target, 230 k events with the T2K replica target and 80 k events without 
target (empty target events) were registered during the 2007 run. Using these data 
interaction cross sections and charged pion spectra in p-C interactions at 31 GeV/c 
were measured^!. Such measurements are required to improve predictions of the 
neutrino flux for the T2K long baseline neutrino oscillation experiment in Japan ^H}. 

For normalization and cross section measurements we adopted the same procedure 
as the one developed by the NA49 Collaboration^!. The measured inelastic and 
production cross sections in p-C interactions at 31 GeV/c are 257.2 ± 1.9 ± 8.9 mb 
and 229.3 ± 1.9 ± 9.0 mb, respectively. 

Crucial for this analysis is the identification of the produced charged pions. De- 
pending on the momentum interval, different approaches have been adopted, which 
lead also to different track selection criteria. The task is facilitated for the negatively 
charged pions, by the observation that more than 90% of primary negatively charged 
particles produced in p-C interactions at this energy are 7r~, and thus the analysis of 
7r~ spectra can also be carried out without additional particle identification. 

Three analysis methods were developed to obtain charged pion spectra. These 
are: 1) analysis of 7r~ mesons via measurements of negatively charged particles (h~ 
analysis); 2) analysis of 7r + and 7r~ mesons identified via dE/dx measurements in 
the TPCs (dE/dx analysis at low momentum) and 3) analysis of ir + and tt~ mesons 
identified via time-of-flight and dE/dx measurements in the ToF-F and TPCs, re- 
spectively (tof — dE/dx analysis). Each analysis yields fully corrected pion spectra 
with independently calculated statistical and systematic errors. The spectra were 
compared in overlapping phase-space domains to check their consistency. Comple- 
mentary domains were combined to reach maximum acceptance. 

The agreement between spectra obtained by different methods is, in general, better 





Figure 5: Laboratory momentum distributions of 7r + mesons produced in production p-C interac- 
tions at 31 GeV/c in different intervals of polar angle (9). The spectra are normalized to the mean 
tt + multiplicity in all production p-C interactions. Error bars indicate statistical and systematic un- 
certainties added in quadrature. The overall uncertainty (2.3%) due to the normalization procedure 
is not shown. Predictions of hadron production models, FLUKA2008 (solid line), URQMD1.3.1 
(dashed line) and VENUS4.12 (dotted line) are also indicated. 



than 10%. In order to obtain the final spectra consisting of statistically uncorrelated 
points the measurement with the smallest total error was selected. 

Inclusive production cross sections for negatively and positively charged pions 
are presented as a function of laboratory momentum in 10 intervals of the laboratory 
polar angle covering the range from up to 420 mrad^l. The final spectra for tt + are 
plotted in Fig. [5j For the purpose of a comparison of the data with model predictions 
the spectra were normalized to the mean 7r + multiplicity in all production interactions. 
This avoids uncertainties due to the different treatment of quasi-elastic interactions 
in models as well as problems due to the absence of predictions for inclusive cross 
sections. 

As a first application of these measurements, it is interesting to compare the 7r + 
spectra in p-C interactions at 31 GeV/c to the predictions of event generators of 
hadronic interactions. Models that have been frequently used for the interpretation 



of cosmic ray data, i.e. VENUS4.12 FLUKA2008 ^ an d URQMD1.3.1 46) were 
selected. They are part of the CORSIKA framework for the simulation of air 
showers and are typically used to generate hadron-air interactions at energies below 
80 GeV. In order to assure that all relevant settings of the generators are identical to 
the ones used in air shower simulations, p-C interactions at 31 GeV/c were simulated 
within CORSIKA in the so-called interaction test mode. 

The results are presented in Fig. [5j The URQMD1.3.1 model qualitatively dis- 
agrees with the data only at low momenta (p < 3 GeV/c) and polar angles below 
about 140 mrad. The VENUS4.12 and FLUKA2008 models follow the data trend in 
all measured polar angle intervals. 

The data presented here and in Ref. already provide important information 
used for improved calculation of the T2K neutrino flux. Meanwhile, a much larger 
data set with both the thin (4% Ai) and the T2K replica carbon targets was recorded in 
2009 (about 6 M triggers with the thin target and 2 M triggers with the replica target) 
and 2010 (about 10 M triggers with the replica target) and is presently being analyzed. 
This will lead to results of higher precision for pions and extend the measurements 
to other hadron species such as charged kaons, protons, Kg and A. The new data 
will allow a further significant reduction of the uncertainties in the prediction of the 
neutrino flux in the T2K experiment. 

4. Conclusions 

The HARP experiment has already made important contributions to the cross- 
section measurements relevant for neutrino experiments. 

The MIPP experiment has nearly finalized measurements of hadron production 
from the NuMI target used in the MINOS experiment. 

First measurements of charged pion production in p-C interactions at 31 GeV/c 
released recently by NA61 (SHINE) are of significant importance for a precise pre- 
diction of the J-PARC neutrino beam used for the first stage of the T2K experiment. 

All three experiments provide also a large amount of input for validation and 
tuning of hadron production models in Monte-Carlo generators. 

Clearly, hardon production studies is a must for precision neutrino experiments. 
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